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Mean-field study of charge order with long periodicity in 0-(BEDT-TTF)2X 

Masato Kaneko * and Masao Ogata 
Department of Physics, University of Tokyo, Tokyo 113-0033 

Charge ordering phenomenon in ^-(BEDT-TTF)2X is studied by using 1/4- filled extended 
Hubbard model on an anisotropic triangular lattice through mean-field approximation. It is 
found that a metallic charge ordered state with 3-periodicity on lattices (3-fold type charge 
order) is realized in the realistic parameter region where the nearest neighbor Coulomb inter- 
action Vij is nearly isotropic. This 3-fold state survives up to high temperatures by estimating 
the free energy. Insulating state with diagonal or vertical-type charge ordering appears as 
increasing anisotropy of Vij. Considering the effect of the structural distortion observed in 
S-(BEDT-TTF)2RbZn(SCN)4 at low temperatures, horizontal-type charge ordered insulating 
phase tends to be stabilized and the region where 3-fold type exists becomes narrower. Our re- 
sults are consistent with the metallic state with long periodic charge order which can be related 
to 3-fold type in 6l-(BEDT-TTF)2RbZn(SCN)4 at high temperatures and insulating state with 
horizontal charge order at low temperatures. For 5-(BEDT-TTF)2CsZn(SCN)4, we speculate 
several kinds of charge-ordered states are energetically competing with each other leading to 
an inhomogeneous state at low temperatures. 

KEYWORDS: charge order, BEDT-TTF, organic conductor, extended Hubbard model, triangular lattice, 
Hartree approximation 



1. Introduction 

The family of organic conductors, (BEDT-TTF)2X (in 
short (ET)2X), which consist of two-dimensional stack- 
ing layers of donor BEDT-TTF (ET) molecules has a 
variety of electronic properties, such as metal-insulator 
transition, superconductivity and charge ordering, de- 
pending on temperature, pressure and the details of its 
anion molecule X.^'^ This indicates that the effect of 
strong correlation is important and various states are 
competing with each other. In particular, the value of 
the nearest neighbor Coulomb interaction is considered 
to be relatively large in organic compounds and crucial 
for charge ordering. The various properties observed ex- 
perimentally have been discussed on the basis of theo- 
retical models which regard the ET molecule as a single 
site.^'^"^ 

In this paper, we focus on the 0-type and ^^-type 
(ET)2X which have interesting behaviors mentioned 
later. Here 6 and 9d indicate specific lattice structures 
leading to different transfer integrals and intermolecu- 
lar interactions. It is considered that the charge ordering 
phenomena seen in these materials can be theoretically 
explained in the 1/4-filled extended Hubbard model in- 
cluding nearest neighbor Coulomb repulsion V between 
ET molecules.^ For the 0-type (ET)2X, the model con- 
sists of an anisotropic triangular lattice described in 
Fig.l. Two kinds of nearest neighbor Coulomb repulsions 
and transfer integrals between molecules are denoted by 
the subscript p and c corresponding to diagonal and ver- 
tical direction, respectively. This system is characterized 
by the dihedral angle cjP' between molecules which de- 
pends on pressure or anion X and controls anisotropy 
expressed as the ratio of nearest neighbor Coulomb re- 
pulsion Vp/Vc and transfer integral tp/tc- According to 




Fig. 1. The structure of 9-type (ET)2X. The two types of transfer 
integrals and nearest neighbor Coulomb interaction are denoted 
by p and c for diagonal and vertical direction, respectively. The 
dihedral angle </> between molecules is also shown. 



the angle, various physical properties appear. 

For example, 6'-(BEDT-TTF)2RbZn(SCN)4 shows a 
metal insulator transition at Tmi = 200K under slow 
cooling speed. This transition is accompanied by a struc- 
tural distortion which gives 6d structure below Tmi- Fur- 
thermore, a horizontal stripe charge order is observed. 
In the lower temperature (T < 25K), a spin gap behav- 
ior is observed. These properties were clarified based on 
various experimental measurements.^' ^"^^ However, a re- 
cent X-ray scattering experiment^ has shown that, even 
in the metallic phase, some kinds of short-range charge 
order with longer periodicity exists. This charge order- 
ing is different from the one seen in the insulating phase 
below Tmi (i-e., ^^-phase). Moreover, the NMR stud- 
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cates a feature of insulators. In contrast, 0-(BEDT- 
TTF)2CsZn(SCN)4 shows metallic resistivity down to 
20K, below which resistivity exhibits sudden increase. 
A sign of growing horizontal charge order is observed 
in X-ray measurements.^^ Furthermore, large dielectric 
constant and giant nonlinear conduction are observed at 
low temperatures, suggesting an appearance of inhomo- 
geneous charge ordered state. This is also consistent 
with NMR and EPR results^-^'^^ and can be related to 
the' optical conductivity spectra showing a sudden peak 
at low frequency as lowering temperature.^'^ 

Having these experiments in mind, we study the 1 /4- 
filled extended Hubbard model in a mean-field approx- 
imation. We find that a charge ordered state with long 
periodicity is realized in ^-phase when the frustration in 
the nearest-neighbor Coulomb repulsions increases. We 
call this state as a 3-fold state. Furthermore, we show 
that the various charge ordered states have similar en- 
ergies, which will be the origin of the appearance of in- 
homogeneity observed experimentally. In the ^^-phase, 
we show the competition between the 3-fold state and a 
horizontal stripe state. 

The extended Hubbard model has been studied in var- 
ious contexts. On a square lattice where Vc and tc are 
equal to zero, it has been confirmed that a charge or- 
dred state is stabilized above critical values of V. 24-26 
However, the actual ratio of Vp/Vc in 0-(ET)2X com- 
pounds has been estimated to be comparable, namely 
Vp/Vc — Thus, it is necessary to investigate an 
anisotropic triangular lattice where a strong frustration 
effect due to the lattice structure is expected with re- 
spect to the charge degree of freedom.^ Using realistic 
transfer integrals, Seo found various patterns of charge 
order in the ground state by using a mean-field approx- 
imation.^ However, he used rather small unit cells and 
did not take into account of states with longer period- 
icity. Recently, Mori discussed several states with long 
periodicity in the atomic limit, i.e., the kinetic energy is 
completely neglected.^ In this paper, we study the unre- 
stricted Hartree approximation in the extended Hubbard 
model with kinetic energy using the realistic parameters 
for ^^-type and ^^^-type (ET)2X compounds. This enables 
us to study both metallic and insulating charge ordered 
state with long periodicity. 

In the next section, the model and mean-field approx- 
imation is introduced in order to discuss the charge or- 
dering phenomena. The results of calculations which con- 
tain various kinds of charge ordered states are presented 
in § 3. The 6'-phase and 6d-pha.sc, realized in 6'-(BEDT- 
TTF)2RbZn(SCN)4 at high and low temperatures, are 
studied in § 3.1 and 3.2, respectively. In § 4, we discuss 
the obtained results in view of the relation to the exper- 
imental facts for the actual materials. 

2. Method and Formulations 

We study the quarter filled extended Hubbard model 
on an anisotropic triangular lattice, 



<ij>(T i <ij> 



where U and Vij are the on-site and the nearest neigh- 
bor Coulomb repulsion, respectively, which are crucial for 
charge order. Vij and tij depend on the direction between 
the neighbor molecules, as shown in Fig.l where the lat- 
tice structure of the ^-type (ET)2X and two types of Vij 
and tij are shown. We will also study the ^^-type case 
later, which has a little complicated Vij and tij (Fig. 8). 
In eq.(l), < ij > represents the summation over the pairs 
of neighbor sites, a, riia and {cia) denote up or down 
spin, number operator {rii = riii +nii) and the creation 
(annihilation) operator for an electron with spin a at ith 
site, respectively. 

The mean-field approximation at zero temperature of 
the above model was carried out by Seo,'' by decoupling 
the interaction terms as Uii^nn (^it)^ii + '^iT(^ii)- 
Although various patterns of charge order could be ob- 
tained, Seo studied rather restricted charge order pat- 
terns. In order to include the possibility for charge or- 
der with long periodicity implied by recent experiments, 
we carry out an unrestricted Hartree method using the 
larger number of sites in a unit cell up to 36 sites both 
at zero and at finite temperatures. 

Through mean-field approximation, eq.(l) can be writ- 
ten in k space as follows, 



kcr Im 



(2) 

where /' and m are the index of sites in a unit cell. By 
diagonalizing the matrix (2), the self-consistency equa- 
tions for the electron density at i th site, [uir,) are solved. 

To determine which pattern of charge order is the most 
stable among the various solutions, the ground state en- 
ergy E and the free energy F at finite temperatures are 
calculated and compared with each other. These quanti- 
ties per unit cell are given by. 



cell 



<lm> 



Vlm{ni){nrn), 



(3) 



^'cell 



-/Vcell 



ikcr 



-^C/(n;|)(ni|) - ^ Vim{ni){nm), (4) 



<lm> 



(5) 



where £^;k<7 and uimka are the eigenvalue and the element 
of unitary matrix obtained by diagonalization of (2), re- 
spectively. nF{Ei]s_^) is the Fermi distribution function 
■^^j^^niEi^^-,,) and /X is the chemical potential determined 
from (5). Nioi and A'^ccii are the total number of electrons 
and unit cells, respectively. 
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3. Results 

3.1 e-type 

First, let us study the case of 9-type (ET)2X as shown 
in Fig.l. Coulomb interaction U, diagonal hopping tp and 
the ratio Vc/U are fixed at [/ = O.TeV, tp = O.leV and 
Vc/U = 0.3, respectively which are considered as typi- 
cal values for organic conductors. To investigate the 
dependence on the anisotropy, the ratio of the trans- 
fer integral and nearest neighbor Coulomb interaction 
are varied as \tc/tp\ =0—1 and Vp/Vc = — 1, re- 
spectively. This choice of parameters is the same as in 
ref.^ As for realistic values, |tc/ip| is estimated to be 
almost equal to 0.1,^'"^^ while Vp and Vc are compa- 
rable. This means that the transfer integrals are sim- 
ilar to the square lattice, while the nearest neighbor 
Coulomb interactions have triangular lattice character. 
In this model, noninteracting energy dispersion is given 
by -Ek = 2tcC0s{ky) ± 4tp cos(fca;/2) cos(A:y/2). The sign 
of tp is irrelevant, while that of tc is fixed to be positive 
for the quarter-filled hole system. 

We find five charge order patterns as shown in Fig. 2, 
where the charge rich and poor sites as denoted by 
solid and open ellipsis, respectively. We have tried many 
other patterns, but they are not stabilized or have a 
rather higher energy. Fig. 2(a) shows a pattern with a 
unit cell with three lattice sites which was not consid- 
ered by Seo.^ We call this state as "3- fold" pattern in 
the following. The other four states are stripe states 
and we call them as (b)vertical (c)horizontal (d)diagonal 
(e)vertical-l-diagonal, respectively. 

The obtained phase diagram on tc/tp — Vp/Vc plane 
for U = 0.7eV, VJU = 0.3, U/tp = 7 is shown in Fig.3. 
Here, the lowest-energy state among the various states 
is shown. The diagonal-1 and diagonal-2 (and 3-fold-l 
and 3-fold-2) have different spin configurations which we 
discuss shortly. We find that the 3-fold pattern is stabi- 
hzed in the large Vp/Vc region. The vertical stripe state 
obtained in^ has a slightly higher energy than the 3-fold 
pattern in this region. The phase diagram in Fig.3 shows 
that the frustration in Vij leads to the 3-fold pattern. 

The hole density at each site depends on the strength 
of Coulomb interaction and transfer integral. Their val- 
ues in the ground state range from 0.7 to 0.9 at rich sites 
and from nearly to 0.1 at poor sites, except for the 3- 
fold type where they are around 1.0 and 0.25 at rich and 
poor sites, respectively. 

The diagonal-1 and diagonal-2 states are insulating 
states in the parameter region where calculations were 
carried out, since the Fermi energy is located in the 
band gap. On the other hand, the 3-fold type states are 
found to be metallic, because the Fermi energy necessar- 
ily crosses at least one band among the six bands in the 
3-fold type states at quarter-filling. 

Next, let us discuss the spin configurations. As shown 
in Fig. 4, not only the charge degrees of freedom but 
also the spin degrees of freedom are ordered. In the 
stripe states, antiferromagnetic order is obtained along 
the stripe due to the strong on-site Coulomb interac- 
tion leading exchange coupling Jy. However, the relative 
phase of antiferromagnetism between neighboring stripes 



(a) 3fold (b) vertical 




(c) horizontal (d) diagonal 




(e) vertical+diagonal 




Fig. 2. Obtained various patterns of charge order. The charge 
rich and poor sites are denoted by sohd and open eUipses, re- 
spectively. 





' ferro ' / 


1 1 


0.8 




3fold-2 


0.6 


diagonal-2 




0.4 






0.2 


diagonal-1 


/ 3fold-1 





1 1 


/ 1 



0.2 0.4 0.6 0.8 1 

Vp/Vc 

Fig. 3. The phase diagram on tc/tp — Vp/Vc plane for U = 
O.leV, Vc/U = 0.3, U/tp = 7. The differences between diagonal- 
1 and -2 and between 3fold-l and -2 are spin configuration as 
explained in Fig. 4. 

depends on the inter-chain coupling. The diagonal-1 has 
antiferromagnetic spin order along the c axis (i.e., verti- 
cal direction) between stripes, while the diagonal-2 has 
ferromagnetic spin order as shown in Fig. 4(c) and (d). 
These two states have very close energies and this will be 
due to the frustration effect of the triangular lattice. On 
the other hand, in the 3-fold type charge order, the spins 
on the charge-rich sites have the same direction, as shown 
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(c) diagonal-1 



(d) diagonal-2 




Fig. 4. 3-fold and diagonal type charge order appeared in Fig.3. 
The grey and open ellipses represent charge rich and poor sites, 
respectively. The direction of spin at each site is shown by up- 
arrow or downarrow. The difference between the diagonal-1 and 
the diagonal-2 is the spin configuration at charge rich sites. The 
former has antiferromagnetic spin order along c axis between 
stripes, while the latter has ferromagnetic spin order. 3fold-l is 
ferrimagnetic, and 3fold-2 is ferromagnetic. 



in Fig. 4(a) and (b). The 3-fold-l has antiferromagnetic 
spin order between charge rich and poor sites leading to a 

ferrimagnetic state, while the 3-fold-2 has ferromagnetic 
spin order. We find that the 3-fold-l is the most stable in 
the large Vp/Vc and small \tc/tp\ region. This is because 
the gain of the antiferromagnetic exchange energy is the 
largest among the various spin configurations. In the 3- 
fold-1, a charge-rich site is surrounded by six charge-poor 
sites with opposite spin direction. 

The phase diagram in Fig.3 shows that the charge or- 
der pattern is insensitive to the ratio of transfer \tc/tp\, 
although the ordering of spin is varied and ferromagnetic 
spin configuration is favored as tc/tp approaches 1. We 
think that this tendency to the ferromagnetism is due to 
the van Hove singularity. When tc{> 0) becomes larger, 
the van Hove singularity approaches the Fermi energy at 
quarter filling. This leads to a large density of state which 
causes the ferromagnetism. For the triangular lattice, it 
has been shown that the ferromagnetism appears in the 
large U region irrespective of the electron density. 

In Fig. 5, we show the ground-state energy difference 
between the 3-fold and the other stripe-type states for 
U = O.leV, Vc/U = 0.3, U/tp = 7 and \tc/tp\ = 0.1. We 
can see that the 3-fold metallic state is stabilized in the 
region of Vp/Vc > 0.64, while insulating state with di- 
agonal stripe charge order is realized for Vp/Vc < 0.64 
where the anisotropic effect by Coulomb interaction is 
strong. This is reasonable because the ground state with 
charge order along diagonal is energetically more favor- 
able when Vp/Vc is small enough. In the same way, the 
vertical type charge order has a lower energy than the 
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Fig. 5. The energy difference between the 3-fold and the other 
stripe type states for TJ = O.TeV, Vc/C/ = 0.3, ?7/tp = 7 and 
Kc/ip| = 0.1. The energy of the 3-fold type state is choose to be 
zero. 



diagonal type when Vp/Vc becomes larger. However, the 
3-fold type has a much lower energy in this region, as 
shown in Fig. 5. 

The behavior at finite temperatures is also investigated 
by calculating the mean- field free energy F. As seen in 
Fig. 6, we find that the 3-fold state gives the lowest free 
energy in a wide range of temperature when Vp/Vc — 1- 
The other charge ordered states have higher free ener- 
gies which are close to each other. When Vp /Vc is small, 
the free energy of the state with diagonal charge stripe 
becomes the lowest. Interestingly, in the intermediate 
values of V^/K, we find that a first-order phase tran- 
sition between the diagonal stripe state and the 3-fold 
state at a finite temperature. For example, the free en- 
ergies at Vp/Vc = 0.6 are shown in Fig. 7. Apparently 
the energies of various states become close to each other, 
and the diagonal-1 state is replaced by the 3- fold state 
at T/tp = 0.2 as the temperature increases. This phase 
transition shows the higher entropy gain of 3-fold phase 
represented in the gradient of free energy. 

3.2 Od-type 

In this subsection, we discuss the ^d-phase which is 
realized in some of (ET2)X materials at low tempera- 
tures. The crystal structure of the (^,;-phase is not so 
simple as the 0-phase due to the lattice distortion. There- 
fore, different parameters from those of the 0-phase need 
to be introduced for calculations. In the 0af-phase, the 
values of hopping parameter tij have been estimated 
through extended Hiickel method based on experimen- 
tal data measured on the 6'-(BEDT-TTF)RbZn(SCN)4 
at T < Tmi,^'^'^ which are shown in Fig. 8 together with 
its structure. The ratio of Vp /Vc is expected to be almost 
the same as in the O-phasc since the distance between the 
molecules does not change so much. We carry out calcu- 
lations by using the set of parameters Uj given in Fig. 8 as 
varying the ratio Vp/Vc with U = 0.7eV and Vc/U = 0.3. 

The energy difference between 3-fold, horizontal, ver- 
tical and diagonal- type is shown in Fig. 9 as a function of 
Vp/Vc- Here, we plot the lowest energy state for each type 
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phase is stable widely for finite temperatures. 
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Fig. 7. Free energy for Vp/Vc = 0.6. The other parameters are 

fixed at f/ = 0.7eV,Vc/U = 0.3, U/tp = 7 and \tc/tp\ = 0.1. A 
first-order phase transition between the diagonal-type insulating 
state and the 3-fold-type metallic state occurs at T/tp = 0.2. 



of charge order among the obtained states with differ- 
ent spin configurations. It is found that, while the 3-fold 
charge ordered state is stabilized for Vp/Vc — 1, the hori- 
zontal charge ordered state is stabilized for Vp/Vc < 0.92. 
Considering the horizontal state is never stabilized in the 
case of the ^-type, the lattice distortion plays an im- 
portant role for the horizontal stripe charge order. This 
tendency is also found in the calculation by exact diago- 
nalization with on-site electron-phonon effects. 

Among the obtained horizontal- type states, the one 
with antiferromagnetic charge stripe formed along lines 
connected by the large transfer tp4 has the lowest energy. 
This is due to the effect of exchange coupling between the 
nearest neighbor spins (~ tp^/U). Our result is consistent 
with the result by Seo where different values of transfer 
integrals were used. On the other hand, in the 3-fold type 
realized for the large Vp/Vc, there are three patterns of 
charge distribution in a unit cell with one kind of spin 
configuration. They are within almost the same energy 
level and have a tendency with charge rich and poor sites 
having opposite spin, although this is not so simple as in 
Fig. 4(a) and has rather complicated order. 

4. Discussions 

In this section, we discuss the results described in the 
previous section in relation to the charge ordering phe- 




Fig. 8. Od phase is characterized by the six kinds of transfer 
integrals tpi ~ tp4, td and tc2- These values arc estimated 
through extended Hiickel method as follows, tpi = 16.9, tp2 = 
-6.5, tp3 = 2.2, tp4 = -12.3, tci = 1.5 and tc2 = 5.2 (lO-^eV). 
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Fig. 9. The energy difference between 3-fold, horizontal, vertical 
and diagonal-type. The energy of the 3-fold type is chosen to be 
zero. Parameters are fixed at [/ = 0.7eV and Vc/U = 0.3. 



nomena observed experimentally in 0-(BEDT-TTF)2X. 

The 3-fold type solution obtained in the present cal- 
culation can be related to the X-ray scattering exper- 
iment^ in the metallic phase above Tmi in ^-(BEDT- 
TTF)2RbZn(SCN)4. From the calculation of the free en- 
ergy, this 3-fold metallic state is the most stable even at 
high temperatures in the ^-phase. Our calculation shows 
that the parameter Vp/Vc — 1 is important for the stabi- 
lization of the 3-fold state. This value is consistent with 
estimation.^ 

However, with respect to the spin degrees of freedom, 



6 J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



the lowest-energy state in our approximation is ferri- 
magnetic, while paramagnetic susceptibility has been ob- 
served experimentally.^' This will be due to the mean- 
field approximation where the effect of quantum fluctu- 
ation is neglected. Thus, in order to discuss the spin de- 
grees of freedom precisely, more appropriate treatment 
is required, although the charge degrees of freedom are 
comparable with the actual measurements. 

Below the metal-insulator transition, 0-(BEDT- 
TTF)2RbZn(SCN)4 changes its structure from to 0^- 
type. For this 9d case, we find the stabilization of the hori- 
zontal stripe which is consistent with experiments.^' ^^"^^ 
The horizontal stripe charge order is observed not only 
in X=RbM' but also in X=T1M' (M' = Zn, Co) in the 
0d-phase.^^ Our result suggests that, in addition to the 
ratio Vp/Vc, the network of transfer integral tij plays 
an important role for the stabilization of the charge 
order patterns. The magnetic susceptibility by experi- 
ments can be fit for the the Bonner-Fisher curve for the 
one-dimensional Heisenberg model. ^ Although the mean- 
field calculation gives an antiferromagnetic spin order 
along the horizontal stripe, it is natural to have a one- 
dimensional Heisenberg spin fluctuation at finite tem- 
peratures when the quantum fluctuations are taken into 
account. 

In the present paper, we did not discuss the phase 
transition between 6 and ^^-type structure observed in 
^-(ET)2RbZn(SCN)4. For this purpose, the elastic ener- 
gies, or lattice distortion energies should be taken into 
account, which is beyond the scope of the present paper. 
Generally, the symmetry in the ^-phase can be broken to 
a lower-symmetry as the temperature decreases. Thus, it 
is natural to expect Od-phase with insulating horizontal 
stripe at low temperatures. 

Experimentally, the other type of charge ordered 
states, i.e., diagonal and vertical-type have been actu- 
ally found. For example, ^m-TlZn exhibits the former 
type.^^ Here, 9m represents a monoclinic structure in 
contrast to an orthorhombic one in 6 and ^^-phase. On 
the other hand, 6l-(BDT-TTP)2Cu(NCS)2 having simi- 
lar structure to the 6'-typc exhibits a checkerboard type 
charge order. Thus, the implication of various compet- 
ing states makes the system complicated and causes in- 
teresting physical properties. 

In the case of i9-(BEDT-TTF)2CsZn(SCN)4, the inho- 
mogencous emergence of charge ordered insulating do- 
mains is suggested.^" These domains seem to coexist with 
other metallic domains. On the other hand, the X-ray 
diffraction study shows that a horizontal charge order 
grows blow 50K.^^ These experimental observations will 
be related to our result where 3-fold metallic state and 
horizontal insulating charge ordered state are competing 
with each other energetically. In the phase diagram for 9- 
(ET)2X families, 6'-(ET)2CsZn(SCN)4 is located near the 
boimdary between 9 and the ^^-structure.^ Therefore, it 
is possible that the strength of coupling to the lattice 
grows only in some local domains at low temperatures. 
As a result, competing phases can gradually form in each 
domain, although the global structural phase transition 
to ^d-type as in ^-(ET)2RbZn(SCN)4 does not occur. 

In summary, we investigated the charge ordering phe- 



nomenon in 0-(ET)2X by using 1/4-filled extended Hub- 
bard model on an anisotropic triangular lattice through 
mean-field approximation. It is found that metallic state 
with 3-fold type charge order is realized in the 0-phase 
in the parameter region where the nearest neighbor 
Coulomb interaction Vij is nearly isotropic. This 3-fold 
state survives up to high temperatures. Insulating state 
with diagonal or vertical-type charge ordering appears 
as increasing anisotropy of Vij. On the other hand, in 
the 9d-phi\sc. horizontal-type charge ordered insulating 
phase is stabilized and the region where 3-fold type ex- 
ists becomes narrower. 

These results are consistent with the metallic state 
with long periodic charge order which can be related 
to 3-fold type in 6'-(ET)2RbZn(SCN)4 at high temper- 
atures. The insulating state with horizontal charge order 
is enhanced in the ^^^-phase at low temperatures. For 9- 
(ET)2CsZn(SCN)4, this effect is considered to play an 
important role within local domains, leading to the in- 
teresting competition with metallic state at low temper- 
atures. 
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